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DESIGN OF DRUGS INVOLVING 
RECEPTOR-LIGAND-DNA INTERACTIONS 

Cross-reference to Related Cases 

This application is a continuation-in-part of 
application serial number 08/158,689 filed on November 26, 1993 
which is hereby incorporated by reference in its entirety. 

Field of the Invention 

The present invention relates generally to rational 
drug design, in particular design of biologically active molecules 
using pharmacophores prepared according to the present 
invention. 

Background of the Invention 

Why certain chemical structures and not others are 
present in nature has been a recurring question raised by scientists 
since the first organic natural products were characterized. Of 
equal interest has been elucidating what structural features within 
any given class of organic molecules are responsible for 
biological activity. Historically, the lack of satisfactory answers 
to both questions has relegated the development of biologically 
active molecules either to serendipity or to exhaustive synthesis 
and biological testing of large numbers of compounds. This 



SUBSTITUTE SHEET (RULE 26) 



wo 95/14791 PCT/US94/13765 

2 

frustration is particularly evident in the pharmaceutical industry 
where the development of drug agonists and antagonists is often 
time consuming, tedious and expensive. 

This picture is beginning to change as more 
5 informadon is derived from modem molecular modeling 

techniques including characterization of the active sites in 
enzymes arid the hgand binding sites in receptors. Over the past 
15 years, another approach has emerged based upon a series of 
discoveries made with molecular models, wherein biologically 
10 active small molecules have been found to possess complementary 

stereochemical relationships with gene structure. This approach 
was first described in U.S. Patent No. 4,461,619 to Hendry, et al., 
which is incorporated herein by reference. This simple molecular 
modeling technology was developed from observations first 
15 reported in 1977 of structural relationships between small 

molecules and nucleic acids, as described by Hendry, et al., J, 
Steroid Biochem. Molec. Biol 42:659-670 (1992); Copland, et 
al., /. Steroid Biochem. Molec. Biol 46:451-462 (1993); Hendry 
and Mahesh, 7. Steroid Biochem. Molec. Biol 41:647-651 (1992); 
20 Witham and Hendry, /. Theor. Biol 155:55-67 (1992); Hendry 

and Mahesh, /. Steroid Biochem. Molec. Biol 39:133-146 (1991); 
Hendry, /. Steroid Biochem. 31:493-523 (1988); Lehner, et al., 
Molec. Endocrinol (1987); Hendry, et al., /. Steroid 

Biochem. 24:843-852 (1986); Uberoi, et al.. Steroids 45:325-340 
25 (1985); Bransome, et aL, /. Theor, Biol 112:97-108 (1985); 

Hendry, et al., Proc. Natl Acad. Scl USA 78:7440-7444 (1981); 
and Hendry, et al., Perspect. Biol Med. 27:623-651 (1984), all of 
which are incorporated herein by reference. 
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25 



The essential ingredient of all genes is a single, well 
defined polymer, deoxyribonucleic acid (DNA). DNA is a 
remarkably uncomplicated molecule composed of recurring 
sugar-phosphate units attached to one of four possible bases: 
adenine (A), thymine (T), cytosine (C) or guanine (G). The 
simplicity of gene structure is further evident in the Watson and 
Crick base pairing scheme of double-stranded DNA (A with T 
and C with G), and the helical chirality (handedness) dictated by 
the absolute configuration of the sugar D-deoxyribose. Gene 
structure could conceivably be composed of many other chemical 
units, for example, other sugar stereoisomers such as L- 
deoxyribose or sugar homologs related to D-glucose. 

The products of gene structure, proteins, are also 
simple, ubiquitous molecules. Nature limits the structure of 
proteins by constructing them from only twenty basic units, the 
amino acids; protein chirality is constrained by the absolute L- 
configuration of the amino acids. As in the case of nucleic acid 
subunits, a wide range of structural alternatives are possible for 
protein amino acids. Examples include changes in the chiraUty of 
a given amino acid side chain (e.g., D-isoleucine), rearrangements 
in the pattern of atoms (e.g., the t-butyl isomer of isoleucine) or 
the addition of atoms (e.g., pipecoUc acid, a homologue of 
proline). 

Structural constraints are also evident in the 
stereochemistry of low molecular weight natural products. 
Particularly conspicuous are Umitarions imposed by nature on the 
number, size, shape, elemental composition, and chirality of 
biologically active small molecules. For example, the pervasive 
neurotransmitters histamine and serotonin are unique in that 
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alternative structures with changes in the position or composition 
of heteroatoms and/or ring patterns generally do not exist in 
nature. Similarly, many small molecular weight hormones are 
few in number, have recurring structural patterns and possess a 
5 single absolute chirality. 

The source of the pervasive occurrence of 
physicochemical constraints on the structure of naturally 
occurring small molecules lies directly in the structure of the 
proteins which govern both their biosynthesis and bioactivity, i.e., 

10 enzymes and receptors, respectively. Ultimately, however, this 

stereochemical information is contained in the genes. According 
to the basic tenants of molecular biology, the information in DNA 
is replicated with remarkable precision and fidelity into newly 
synthesized DNA. It is also transcribed into RNA and 

15 subsequently translated into protein. 

This scenario, however, presents an apparent 
paradox. While the genetic template ultimately directs which 
proteins and small molecules are synthesized, as well as which 
proteins and small molecules will interact with each other, the 

20 undirectional flow of genetic information during translation 

suggests that DNA structure performs this function without 
recognizing the structure of the small molecule. With few 
exceptions, such as certain antibiotics which bind directly to DNA 
and block transcription, small molecules are not considered to 

25 recognize or interact with the genetic template. Moreover, the 

structures of the molecules that are biosynthesized are thought to 
be unrelated to the structure of the genes. 

In the initial search for structural relationships 
between biologically active natural products and DNA, it became 
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apparent that the two-dimensional structures of DNA base pairs 
were analogous to many classes of small molecules, including 
gibberellic acid, a phytohormone; benzo [a] pyrene oxide, a 
carcinogen; the prostaglandin PGE2; morphine, a narcotic; 
estradiol, a hormone; riboflavin, vitamin Bi2; serotonin, a 
neurotransmitter; and actinomycin, an antibiotic. In addition to 
similarities in size and shape, numerous small molecules contained 
donor/acceptor functional groups at locations where hydrogen 
bonds occurred between the base pairs. When overlaid on. the 

base nair.^ Qrvmf* r*r»mr*r*HT»rlc cii/^Vi on t-Ur* 1 . 

^ , « w**-.*, w u«LA vi.w> , u u wAi C40 uiw piOii 1 ilKJL iliUliC 

gibberellic acid, the steroid hormone estradiol, and 
prostaglandins, contained heteroatoms separated by intemuclear 
distances similar to that of phosphate oxygens on adjacent strands 
of double-stranded DNA. This was particularly evident in 
functional groups attached at the 3 and 17p positions of the 
steroids. 

Using three dimensional Corey-Pauling-Koltun 
(CPK) space filling models, it became apparent that there were 
spaces between base pairs in partially unwound DNA that could 
accommodate a variety of small molecules. For example, 
estradiol could be inserted between base pairs in DNA, and the 
hydroxyl groups at 3 and I7p of estradiol were positioned such 
tiiat tiiey could form hydrogen bonds to phosphate oxygens on 
adjacent strands of DNA. Other steroids, including testosterone 
and progesterone, were also capable of stereochemical insertion 
between base pairs. In each case, complementary donor/acceptor 
hnkages could be formed and the steroid conformed well to the 
topography of the double helix. Attempts to insert any of the 
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non-naturally occurring steroid enantiomers into DNA resulted in 
poor fit in that donor/acceptor linkages were strained or could not 
form, and/or the overall shape of the molecules was incompatible 
with the helical topography of the DNA. 
5 Certain synthetic compounds with hormonal activity 

can also be accommodated within the DNA; in many cases, the fit 
of synthetic compounds such as diethylstilbestrol mimicked that of 
the natural hormone. In addition to mammalian steroids, 
prostaglandins, the insect hormone ecdysone and several 
10 phytohormones were also capable of stereochemical insertion and 

"recognition" by the double helix. In the case of the plant 
hormone gibberellic acid, four stereospecific hydrogen bonds 
could be formed to donor/acceptor positions on the DNA. As 
with the steroids, only the naturally occurring enantiomer of 
15 gibberellic acid conformed to the topography of the double helix. 

One conclusion drawn from these studies is that 
certain chemical shapes, coupled with heteroatom positioning 
compatible with tiiat of the phosphate backbone of DNA and 
hydrogen bond positions of the base-pair template, potentiate 
20 partial or complete recognition between biologically active 

molecules and DNA. 

While it was possible to form complexes between 
DNA and a variety of molecules, amino acids did not initially 
show any clear accommodation to the space between base pairs. 
25 Certain compounds derived from amino acids, for example, 

neurotransmitters, fit into related sites. 

These relationships have been described as a 
stereochemical logic associated with gene structure. The 
stereochemical logic is defined as those unique features of nucleic 
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acid structure which ultimately dictate constraints on molecular 
structure, function, metabolism, and biologic activity. 

The use of molecular modeling as a tool to study 
organic structure has dramatically increased due to the advent of 
5 computer graphics. Not only is it possible to view molecules on 

computer screens in three dimensions but it is also feasible to 
examine the interactions of ligands with various macromolecules 
such as enzymes and receptors, as reviewed by Borman, Chem, 
Eng. News 70:18-26 (1992). An almost baffling array of 
10 software and hardware is now available and vinually all major 

pharmaceutical companies have computer modeling groups which 
are devoted to drug design. 

Modem methods of drug design include studies 
which focus on the binding of a molecule to a protein such as a 
15 polypeptide ligand for a receptor, or a steroid such as an estrogen 

or progesterone for a receptor. Similarly, drugs can be designed 
based upon the interaction of substrates with various enzymes. 
For the most part, however, binding sites in proteins have been 
difficult to characterize. There are many situations where other 
20 mechanisms must be involved to explain the feedback between 

protein regulation and regulation of gene expression. 

What is needed is a method for accurately predicting 
the biological activity of a given compound. The method should 
be easy to perform and should be able to predict both agonist and 
25 antagonist activity. 

Summary of the Invention 

The present invention is a method for identifying 
biological activity of molecules using pharmacophores. 
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According to the present invention, molecules are screened by 
determining the degree of "fit" in the pharmacophore. 

The method according to the present invention can be 
used to identify drugs having increased biological activity or 
5 which have usefulness as antagonists or agonists, including, for 

example, estrogens and anti-estrogens. This method can also be 
used for the following: to predict the fit of compounds into 
nucleic acids, especially DNA; to predict the bioactivity of 
compounds, to screen compounds for toxicity; to design chemical 
10 groups to add to specific sites on molecules to facilitate 

metabolism or render the drug an agonist or antagonist; and to 
create molecules that mimic the activity of the DNA binding 
regions of receptors. 

The present invention also includes pharmacophores 
15 and the method of producing the pharmacophores and the use of 

the pharmacophores in predicting biological activity of a given 
compound. The present invention also includes the design of 
biologically active molecules using the pharmacophore. 

It is therefore an object of the present invention to 
20 provide a method which can be used to design a biologically 

active molecule. 

It is another object of the present invention to 
provide a method to screen and/or evaluate existing compounds 
for toxicological activity. 
25 Still another object of the present invention is to 

provide a method to predict the toxicity of compounds. 

Another object of the present invention is to provide 
a method to predict the toxicity of compounds for specific organs, 
tissues, and ceils. 
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Yet another object of the present invention is to 
provide a method to design compounds that will have particular 
types of biological activities, including, but not limited to, 
hormonal, neurotransmitter, metabolic, genetic, immunologic, 
5 pathologic, toxic, and anti-mitotic activities. 

Still another object of the present invention is to 
provide a method to predict the bioactivity of compounds 
including, but not limited to estrogenic, anti-estrogenic, 
androgenic, anti-androgenic, progestational, anti-progestational, 
^0 mineralocorticoid, retinoid, vitamin D like, thyroid, and 

glucocorticoid bioactivities. 

Yet another object of the present invention is to 
provide a method to create pharmacophores that can be used to 
design compounds such as drugs, hormones, neurotransmitters, 
15 agonists and antagonists more efficiently and economically. 

It is another object of the present invention to 
provide a method to create receptor pharmacophores that are 
molecular models of the portions of receptor molecules that bind 
to nucleic acids. 

^0 Yet another object of the present invention is to 

provide a receptor pharmacophore that can be used to design 
molecules that bind to nucleic acids with different affinity than the 
receptor. 

Yet another object of the present invention is to 
15 provide a pharmacophore that represents the three dimensional 

arrangement of solvent molecules around the ligand 
pharmacophore that binds to nucleic acids. 

Another object of the present invention is to provide 
a pharmacophore that represents the three dimensional 

SUBSTITUTE SHEET (RULE 26) 



wo 95/14791 



PCTAJS94/13765 



arrangement of solvent molecules around the receptor 
pharmacophore that binds to nucleic acids. 

It . is yet another object of the present invention to 
provide a method to create a pharmacophore that is a three 
5 dimensional model of the nucleic acid binding don^iain of the 

receptor and of the ligand molecule that binds to the receptor and 
interacts with the nucleic acid at a different site. 

Another object of the present invention is to provide 
a pharmacophore that represents the three dimensional 
10 arrangement of molecules that can be attached to other 

pharmacophores to modify their biological activity. 

Still another object of the present invention is to 
provide a pharmacophore that represents the three dimensional 
arrangement of molecules that can be attached to other 
15 pharmacophores in order to design sites for enzymatic cleavage. 

These and other objects, features and advantages of 
the present invention will become apparent after a review of the 
following detailed description of the disclosed embodiments and 
the appended claims. 
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Brief Description of the Figures 

Figure 1 is a schematic of a cavity in DNA and the 
numbered sites which accommodate steroid/thyroid/vitanfiin A and 
D ligands: testosterone (1,7); progesterone (2,7); aldosterone 
(2,5,7 and 9); Cortisol (2,3,4,5, and 7); estradiol (1 and 6); 
triiodothyronine (T3) (1 and 6); retinoic acid (6); and 1,25- (0H)2 
vitamin D3 (1 and 6); these are divided into two groups based on 
their interaction with either site 6 or 7. 

Figure 2 is a computer generated space filling stereo 
view of the DNA cavity (A), which fits active estrogens oriented 
by energy calculations into the DNA cavity (B), whereas (C) 
shows the combined active surface of estrogens removed from the 
cavity in DNA that is used to construct the pharmacophore. 

Figure 3 demonstrates a volume contour map 
(yellow) in stereo with dummy atoms (magenta) surrounding the 
active molecules which were used in the construction of the 
pharmacophore (A); the empty pharmacophore (B); fit of the 
highly active estrogen 3,llp,17p-trihydroxy-7a-methylestra- 
l,3,5(10)-triene 11-nitrate ester (hereinafter 7a-methylestradiol- 
liP-nitrate ester reported in Peters et al., 7. Med. Chem. 
32:2306-2310 (1989)) which is acconmiodated completely within 
the pharmacophore (C); and poor fit of the inactive estrogen 9P- 
estradiol which extends appreciably beyond the surface of the 
pharmacophore (D). 

Figure 4 shows a stereo view of volume maps (green 
mesh) and dummy atoms (magenta) for pharmacophores for 
estrogen (A), androgen (B), thyroid (C), anti-estrogen (D), and 
toxicity (E). Dummy atoms are not presented with the toxicity 
pharmacophore volume map. 
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Figure 5, left column, shows the volume maps (green 
mesh) and dummy atoms (magenta) for pharmacophores for 
estrogen (A), androgen (B), thyroid (C), and-estrogen (D), and 
toxicity (E). The right column presents an orthogonal view (90 
5 degree rotation) of the corresponding images in the left column. 

Figure 6 demonstrates the relationship between the 
degree of fit of various compounds to the estrogen 
pharmacophore and the relative uterotropic (estrogenic) activity. 

Figure 7 shows the chemical structures of three 
10 separate molecules, SGI 100, SGI 101, and SGI 102 designed with 

pharmacophore technology. 

Appendix 1 is a QIC computer tape containing data 
files showing volume maps for each of the pharmacophores 
(estrogen, androgen, thyroid, antiestrogen, and toxicity) presented 
15 in Figures 4 and 5; the molecular database containing the dummy 

atoms for each pharmacophore designated by name; and a file 
containing the four color prints (Figures 2-5) submitted with this 
application. 

20 Detailed Description of the Invention 

Many natural products have structures that exhibit 
stereochemical complementarity with nucleic acids, including 
amino acids, phytohormones, cycUc nucleotides, prostaglandins, 
insect hormones, steroid hormones, neurotransmitters, sugars, 
25 peptide hormones, thyroid hormones, pheromones, and vitamins. 

A striking example is the cyclopentanophenanthrene motif 
repeated in all classes of mammalian steroid hormones, such as 
estrogen and progesterone. Another example is the kaurene 



SUBSTITUTE SHEET (RULE 26) 



wo 95/14791 



PCTAJS94/13765 



nucleus containing eight chiral centers which is evident in the 
gibbereUin class of plant hormones. 

There are many ways to examine the stereochemistry 
of various configurations and conformations of nucleic acids. For 
5 example. Silastic polymer models can be constmcted based upon 

computer derived space filling x-ray coordinates to reflect the 
stereochemistry of partially unwound DNA/RNA complexes, 
RNA-RNA complexes, ent-DNA (mirror image of DNA made 
with L-deoxyribose), and apurinic/apyrimidinic sites in DNA. 

10 These cavities reveal a sequence specificity in the fit of many 

molecules. The apurinic/apyrimidinic sites accommodate amino 
acids according to the known genetic code. The plant hormone 
gibberellic acid fits best into the partially unwound site 5'-dTdA- 
3*, 5'-dTdA-3*; members of the mammalian steroid/thyroid 

15 hormone superfamily fit best into 5'dTdG-3*, 5'-dCdA-3\ Each 

class of mamm£dian hormone forms unique stereospecific donor- 
acceptor hydrogen bonds with DNA. The capacity to fit within 
these cavities in the manner of the index biologically active 
molecule correlates with the degree of biologic activity. It is not 

20 possible to fit chiral naturally occurring molecules into ent-DNA. 

However, ent-DNA accommodates the biologic unknown chiral 
enantiomers, such as ent-progesterone. 

According to the present invention, computer 
modeling is used to examine the relationships between compounds 

25 and their fit in helical DNA. Although described herein with 

reference to double stranded, helical DNA, many of the same 
principles are applicable to double stranded RNA, and/or to RNA- 
DNA hybrids. Unless otherwise specified, double stranded RNA 
and DNA are to be considered equivalents as used herein. 
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Computer modeling can be used to view the interactions of 
molecules as well as to measure the energy of a given interaction. 
While a variety of software packages are available for computer 
modeling of molecules, a preferred software package is Sybyl 
software (version 6.03; Tripos Associates, St. Louis, MO) for 
measuring the docking of various small molecular weight ligands 
into DNA. In the examples described herein, the software is run 
on a Silicon Graphics Indigo Extreme equipped with hardware 
stereo, i.e.. Crystal Eyes (StereoGraphics, San Rafael, CA). 
Structures of small molecules are obtained via: the Cambridge 
Crystallographic Database, Lanfield Rd., Cambridge. England; 
construction with the Concord program or from fragment 
libraries and/or modifications of existing x-ray structures 
followed by energy minimization. All energy calculations are 
made using the Sybyl force field with a 1.2 A van der Waals 
parameter for hydrogen, as described by Hendry, et al., /. Steroid 
Biochem, Molec. Biol 42:659-670 (1992) and Hendry, et al, /. 
Steroid Biochem. Molec. Biol. 39:133-146 (1991). Charges are 
calculated using the Gasteiger-Huckel method to include a and K 
bonding. Partially imwound DNA cavities of various double 
stranded dinucleotide sequences are constructed from the Watson 
and Crick canonical B form of DNA by twisting each of the 
fourteen torsional angles along the sugar-phosphate backbone 
while maintaining the structural integrity of hydrogen bonds 
between the base pairs. The y-endo deoxyribose conformation of 
the sugars is employed and Kollman charges are " calculated. 
Within tiiese. constraints, attempts are made to adjust the degree of 
unwinding and the resulting degree of separation of the base pairs 
to accommodate the width of various ligands 
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Each ligand is inserted into the cavity in DNA using 
van der Waals dot, mesh, and space filling surfaces in stereo to 
guide the docking procedure and nunimize steric strain. The 
oxygens of the phosphate groups are permitted to act either as 

5 proton donors or acceptors and oriented to maximize the fit of 

any given ligand. The docking of the molecules is repeated 
several times. The distances between heteroatoms are monitored 
interactively to optimize the direction and distances of potential 
hydrogen bonds. While evaluating fit of compounds in a given 

10 hormone class, attempts are made to insert all candidate ligands 

into the DNA sequence with donor/acceptor linkages conraion to 
the hormone. Donor/acceptor relationships are further 
maximized by appropriate orientation of functional groups on the 
ligands, e.g.. by adjusting the conformation of each structure to 

15 best mimic the fit of the hormone. 

Van der Waals interactions of the candidate ligands 
are optimized with the DNA surface. The force field is used to 
calculate the relative fit of each ligand by assessing the optimum 
favorable energy change resulting from docking the ligand Each 

20 ligand is docked into the DNA and the steric fit is calculated from 

the change in van der Waals energy; the hydrogen bonding fit is 
calculated from the change in electrostatic energy using charges 
on donor hydrogens and acceptor heteroatoms. The energy 
changes are assessed for each ligand. The greater the negative 

25 energy change upon insertion of a given ligand into DNA, the 

more favorable the fit and the more stable the complex. Optimal 
docking is determined when no fiuther increase in fit is observed. 
The most favorable change in energy is selected to compare the 
relative fit of each molecule. The overall fit of each ligand is 
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assessed by adding the change in kcal of the van der Waals and 
electrostatic energies and normalizing the fit to that of the parent 
hormone. 



5 Complementarity of Biologically Active Structures 

Although not wanting to bound by the following 
hypothesis, it is believed that the mechanism of action of 
hormones in the steroid/thyroid superfamily where ligand bound 
receptor is known to regulate hormone responsive genes is 

in Tn#aHi!^f#=»H Kv nbvcir»5»llv tncf*rrina cmfill TnnlpniilftS into DNA, The 

J J — ^ J 55 — - - 

prior art teaches that the ligand causes a specific conformational 
change in the receptor protein which in turn contacts the DNA 
resulting in gene regulation. In contrast, according to the present 
invention, there is a receptor-mediated insertion of the ligand into 

15 DNA. This is consistent with the widely reported lack of 

correlation between hormonal activity and receptor binding for 
estrogenic steroids. According to Brooks, et al.. Recent Advances 
in Steroid Hormone Action, Moudgil, V.K. (ed) 443-466 
(Walther de Gruyter, NY 1987), who made an extensive smdy of 

20 estrogen structure-function relationships, receptor binding is 

essential for target-cell responses elicited oy the steroids. 
However, the affinity of altered estrogens is not direcriy related to 
the character or extent of the response. At the same time, the 
binding of estrogens and other steroids to DNA in the absence of 

25 receptor has been observed to be weak by numerous investigators. 

These include studies of the flat phytoestrogen coumestrol which 
might be expected to easily slide between base pairs. In short, the 
binding of the steroid independentiy to either the receptor or the 
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DNA does not appear to be sufficient to fully explain hormonal 
response. 

Several pieces of evidence indicate that certain 
molecules, such as steroids, that have been shown to fit 

5 remarkably well between base pairs in DNA may elicit strong 

biological responses albeit through relatively weak interactions 
with DNA. Support for this concept is based on in vitro and in 
vivo experimental data and the energy calculations demonstrating 
a correlation between fit of estrogens into DNA and uterotropic 

10 activity, for example, studies of a new class of antitumor agents 

which have led to the discovery of dmgs that are potent but act by 
binding weakly to DNA via intercalation, reported by Lee, et al., 
J, Med. Chem, 35:258-266 (1992). Additional support that 
binding to receptors alone is insufficient to explain activity is 

15 provided by the observation of poor binding of a potent estrogen 

analog, 1 1 P-acetoxyestradiol, which is considerably more active 
than estradiol. 

In general, given that degree of hormonal activity 
does not always correlate well with receptor binding but does 

20 correlate with fit into DNA, it follows that the mode of action of 

steroids involves a stage(s) in which the ligand is recognized by 
both the receptor and the DNA. This conclusion has also been 
reached in studies of anti-androgens. 

It is contemplated as part of the present invention that 

25 binding of the steroid to its receptor serves as a means to 

recognize the general class of hormone (e.g., estrogen versus 
androgen) whereas the type and degree of the fit of the steroid in 
the DNA is largely responsible for governing the magnitude of 
the biologic response. The steroid alone is incapable of proper 



SUBSTITUTE SHEET (RULE 26) 



wo 95/14791 PCT/US94/13765 



18 



insertion into DNA either in vivo or in vitro without the receptor 
and thus would be unable to generate a full hormonal response 
without the receptor. A potentially critical role for the receptor 
upon binding to DNA, possibly in concert with other transcription 
5 factors, is to regulate the physicochemical properties of the site in 

DNA to permit insertion of the steroid, including the degree of 
unwinding, the capacity of heteroatoms to act as either donors or 
acceptors, and the pattern and degree of solvation. This is further 
supported by the decrease in the surface hydrophobicity of the 

^ .1.^ — J^^^ n<ny:] r> fj* £1^ t* try* r> rr ckT* .f or OC 

phosphorylation which enhances binding of the estrogen receptor 
to specific DNA sequences, reported by Denton, et al., 7. BioL 
Chem, 267:7263-7268 (1992), 

The process of a receptor-mediated insertion of 

15 steroids into DNA presumably involves several steps. For 

example, in the case of estrogens, die initial contact of die DNA 
by the steroid-receptor complex could involve a partially exposed 
D ring of the steroid with the A ring still attached to the receptor. 
Analysis of receptor binding data supports this possibiUty. 

20 Stereospecific recognition of die DNA by partial insertion and 

hydrogen bonding of the 17p-hydroxyl of the steroid D ring with 
the 5'-dTdG-3' strand could be followed by complete insertion 
and the recognition within the double helix manifest by die 
linkage of both stereospecific hydrogen bonds. The weak binding 

25 observed for the steroids widi DNA suggests that die steroid/DNA 

complex might be short lived and rapidly reversible. Certain 
estrogen antagonists which are accommodated differently in DNA 
may form longer lived interactions within the site. 
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There are numerous possible scenarios and 
implications of receptor mediated binding of ligands to nucleic 
acids. For example, the lack of effectiveness of certain anti- 
estrogens in tissues which lack estrogen receptor might be due in 
part by the inabiUty of the antagonist to be transported to the 
DNA obviating the insertion of the ligand. Mutation of the site 
which accommodates the estrogens would also result in improper 
recognition of the Ugand and would be predicted to no longer 
permit the gene to be properly regulated either by agonists or 

A AIWA. W AAAJk^AAk t/ W AAAM«.«-ft^«v «wtw*v ,.--,-w" ** q..-^.- 

compound might act, as suggested by the observation that cavities 
in double stranded RNA and RNA-DNA hybrids can 
accommodate various Ugands, for example, estradiol in 5'-rUrG- 
3',5'-dCdA-3\ 



Molecular Modeling 

Molecular modeling was performed with 
Sybyl/Mendyl 5.4 (Tripos Associates, St. Louis, MO) using an 
Evans and Sutherland PS390 graphics computer equipped with a 

20 stereographic viewer. Structures of piperidinedione ligands were 

provided via constmction with the Concord program or from 
fragment libraries followed by energy minimization. Energy 
calculations were made with Sybyl/Mendyl force field and a 1 .2A 
van der Waals parameter for hydrogen. Charges were calculated 

25 using the Gasteigner-Huckel method which includes a and t: 

bonding. Partially unwound DNA was constmcted with 3'-endo 
deoxyribose conformation and KoUman charges. This method has 
been described by Hendry, et al. J. Steroid Biochem. Molec. Biol. 
42:659-670 (1992); Hendry, et al., /. Steroid Biochem, Molec. 
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BioL 39:133-146 (1991), and Hendry et al., /. Steroid Biochem. 
Molec, BioL 49: No. 4-6, pp. 269-280 (1994) the teachings of 
which are hereby incorporated by reference in their entirety. 

The Ugands were inserted into the cavity in DNA 

5 using van der Waals dot surfaces and the stereoviewer to guide 

the docking procedure and minimize any obvious steric strain. 
The distances between heteroatoms were monitored interactively 
to opdnfiize the direction and distances of potential hydrogen 
bonds. Donor/acceptor relationships were further maximized by 

10 appropriate orientation of funcuonal groups on the ligands, e.g., 

by adjusting the conformation of each structure. Attempts were 
made to optimize van der Waals interactions of the candidate 
ligands with the DNA surfaces. The force field was used to assess 
the relative fit of each ligand by quantitating the optimum 

15 favorable energy change resulting from docking the ligand. 

Steric fit was calculated from the change in van der Waals energy, 
the hydrogen bonding fit was calculated from the change in 
electrostatic energy using charges on donor hydrogens and 
acceptor oxygens. The greater the negative energy change upon 

20 insertion of a given ligand into DNA, ±e more favorable die fit 

and die more stable the complex. Docking was completed when 
no further increase in fit was observed. The most favorable 
change in energy was selected to compare the relative fit of each 
molecule. The overall fit of each ligand was assessed by adding 

25 the change in kcal of the van der Waals and electrostatic energies 

and normalizing the value to that of the best fitting molecule 
(100%). It should be noted that while the energies reported here 
were derived from widely used force field calculations, they were 
not empirically derived. Thus, the absolute values in kcal do not 
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have independent experimental significance. At the same time, 
they are valuable indicators of the relative degree of fit into DNA 
of candidate molecules. 

Previous studies using space filling models indicated 
5 diat 3-phenylacetylamino-2,6 piperidinedione was capable of fully 

inserting between base pairs in DNA and forming a stereospecific 
hydrogen bond between the imino proton of the piperidinedione 
ring and a negatively charged phosphate oxygen on the 
deoxyribose-phosphate backbone. Results employing computer 

demonstrated that this compound had favorable van der Waals 
contacts of approximately -17.7 kcal when inserted into DNA 
with an electrostatic energy of approximately -21.7 kcal resulting 
from the stereospecific hydrogen bond (2.7 A) to phosphate. 

15 Increased fit of the ligand was obtained by substituting a para 

hydroxl group on the phenyl ring; this substitution enabled a 
second hydrogen bond to be formed between the hydroxyl group 
and a phosphate oxygen on the adjacent DNA strand. The 
increase in fit measured by energy calculations due to the second 

20 hydrogen bond (2.64 A) was reflected in an additional -24.6 kcal 

in electrostatic energy. Other substitutions which were made on 
the 3-phenylacetylaniino-2, 6-piperidinedione skeleton did not 
significantly increase fit demonstrated by die normalized energy 
calculations for cenain halogenated analogs. 



25 



Synthesis 

The syndiesis of the unsubstituted derivative, 3-N- 
phenylacetylamino-2, 6-piperidinedione has been briefly 
described by Burzynski, et al., Drugs of the Future 10:103 
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(1985), and was used as a general method for the preparation of 
the desired compounds. Appropriate phenylacetic acids were 
reacted with N-hydroxysuccinimide in the presence of N,N- 
dicyclohexylcarbodiimide (DCC) which gave succinimide esters. 
5 The active esters were stable enough to isolate for physical and 

• spectroscopic characterization although the major portions of the 
esters were used for the next reaction without purification. The 
active esters were reacted with L-glutamine in the presence of 
sodium bicarbonate to obtain die glutamine derivatives. However, 

l\j UUC lU UlC UllllUUiUC'O wi-rnu.JLx*A*5 *^m,^j^^*^ w ^ , 

products were directly used for the next reaction. To prepare 
active esters, the glutamine derivatives were again reacted with N- 
hydroxysuccinimide in the presence of DCC to give the active 
esters which without purification were heated at PS-IOOX to 
15 obtain the desired 2,6-piperidinediones in various yields. During 

the heating process the compounds were racemized. 

Biological Evaluation 

These synthetic derivatives were assessed for 

20 biological potency by measuring their growth inhibitory effects 

on various cell Unes using concentrations of 4 nM based on the 
reported IC50 of 3-phenylacetylamino-2,6 piperidinedione in Nb2 
cells. In YAK lymphoma cells, tiie p-hydroxy compound was the 
most active derivative. This compound, p-hydroxy-3- 

25 phenylacetylamino-2,6-piperidinedione was also the most active 

analog when tested in human leukemia (K652) cells. A dose 
response comparison in K562 cells showed tiiat it was more active 
than the unsubstituted compound over the concentration range 
tested (10-5 to 10'^ M). Prolactin stimulated growth of rat Nb2 
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lymphoma cells was inhibited by each of the compounds with p- 
hydroxy-3-phenylacetyiamino-2,6-piperidinedipne manifesting the 
greatest activity. Compound p-hydroxy-3-phenylacetylamino-2,6- 
piperidinedione was more acdve in Nb2 lymphoma cells than 3- 
phenylacetylamino-2,6 piperidinedione over the range tested (10*^ 
M to 10-3 M). 

Further analysis of growth inhibition of the most 
active analog p-hydroxy-3-phenylacetylamino-2,6-piperidinedione 
compared with the parent compound p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione was performed in MCF-7 
(E-3) human breast cancer cells. Both 3-phenylacetylamino-2,6 
piperidinedione and p-hydroxy-3-phenylacetylamino-2,6- 
piperidinedione inhibited estrogen simulated cell growth. In a 9 
day model, p-hydroxy-3-phenylacetylamino-2,6-piperidinedione 
was more active than 3-phenylacetylamino-2,6 piperidinedione 
with IC50 comparable to tamoxifen (3-phenylacetylamino-2,6 
piperidinedione, 3x10*^ M; p.hydroxy-3-phenylacetylamino-2,6- 
piperidinedione, 7x10*6 M, tamoxifen, lxl0"7 M). The open 
chain hydrolysis product of 3-phenylacetylamino-2,6 
piperidinedione, PAG, did not inhibit cell growth even at high 

concentrations (i.e., lO'^M). 

Computer modeling coupled with energy calculations 
confirm that 3-phenylacetylamino-2,6-piperidinedione is capable 
of inserting between base pairs in partially unwound double 
stranded DNA and forming an energetically favorable complex. 
A hydroxyl group placed in the para position of the phenyl ring 
of 3-phenylacetylamino-2,6-piperidinedione enabled formation of 
a second hydrogen bond thereby linking both DNA strands. This 
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added hydrogen bond resulted in a greater fit in the DNA as 
assessed by energy calculations, i.e., 3-phenylacetylaniino-2,6- 
piperidinedione (61%) versus the p-hydroxy derivative p- 
hydroxy-3-phenylacetylaniino-2,6-piperidinedione (100%). 
5 Various substitutions at the para position as well as ortho and 

meta positions witii fluorine and chlorine did not result in a 
significant increase in fit compared to 3-phenylacetylamino-2,6- 
piperidinedione. 

When the analogs of 3-phenylacetylamino-2,6- 

10 piperidinedione were synthesized and examined for the capacity to 

inhibit cancer cell growth, the p-hydroxy derivative, p-hydroxy- 
^ 3-phenylacetylaniino-2,6-piperidinedione, which was predicted to 
be the most active compound based upon fit into DNA was 
consistently found to be the most potent compound. That the 

15 capacity of the hydroxy group of p-hydroxy-3- 

phenylacetylamino-2,6-piperidinedione to form a second 
hydrogen bond to DNA was responsible for the predicted increase 
in activity is further supported by the lack, of increased potency of 
p-fluoro derivative as well as other halogenated derivatives which 

20 were incapable of forming analogous hydrogen bonds. These 

observations demonstrate that among the compounds exannined, a 
correlation exists between degree of fit into DNA and predicted 
biological potency. These findings also support the contention 
that stereochemical complementarity of small molecules with 

25 nucleic acids can be a powerful tool for designing new dmgs. 

The mode of action of these piperidinediones is still 
not proven. One mode of action might involve insertion into 
DNA as suggested by the computer modeling results and the 
observation that DNA synthesis measured by th3miidine 
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incorporation was significantly inhibited upon treatment with 3- 
phenylacetylaniino-2,6-piperidinedione in Nb2 lymphoma cells. 
No covalent adducts of 3-phenylacetylamino-2,6-piperidinedione 
with DNA were detected and the binding was observed to be weak 

5 and reversible, in comparison to classical intercalating drugs. 

Another possible mode of action is suggested by the observation 
that compound p-hydroxy-3-phenylacetylamino-2,6- 

piperidinedione inhibits estrogen stimulated cell growth in MCF-7 
cells, which is comparable to that of the established anti-estrogen 

10 tamoxiien. In contrast to tamoxifen^ however, neither 3^ 

phenylacetylaniino-2,6-piperidinedione nor p-hydroxy-3- 
phenylacetylamino-2,6-piperidinedione exhibited appreciable 
binding for the estrogen receptor. At the same time, direct 
binding of such anti-estrogens to DNA appears to be weak. Taken 

15 as a whole, these findings support the possibility of a weak 

interaction of the piperidinediones with bodi DNA and the 
estrogen receptor, involving a receptor mediated insertion of the 
ligand into DNA. 

20 Development of Pharmacophores * 

Molecular modeling, as described above, facilitates 
the establishment of the best fit of molecules into nucleic acids 
such as double-stranded DNA based on steric and electrostatic 
considerations. Individual molecules, such as estradiol, fit 

25 optimally into specific sites on DNA based -on the location of 

specific nucleotides and the bonding characteristics of individual 
heteroatoms (see Example 2). Molecules that are related to a 
specific molecule such as estradiol but display chemical 
differences will fit into the estradiol site with different degrees of 
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precision: some may fit better and give rise to estradiol agonisdc 
responses while diose widi poor fit display weak estrogenic 
activity. These different molecules may be aligned relative to the 
docking of heteroatoms with heteroatoms on the DNA to optimize 

5 electrostatic interactions. In the creation of pharmacophores 

described below, molecules with activity equal to or greater than 
that of the hormone are chosen for alignment. To date, such 
molecules fit equally well or better than the hormone into DNA 
using the energy calculation methodology described above. 

10 Molecules which do not fit as well into DNA as the parent 

hormone are excluded from inclusion in the construction of the 
phamiacophore. The alignment of the combined surfaces of the 
molecules occupies a specific volume of space thereby forming a 
three dimensional shape. 

15 Pharmacophores are three dimensional arrangements 

of chemical groups related to a given biological activity which 
enables meaningful comparison of molecules exhibiting the same 
biological function (Namto et al., Eur, J. Med. Chem. 20:529-532 
(1985)). Pharmacophores can be derived by simple overlap of 

20 active stmctures or common functional groups in the molecules. 

Without a way to orient the molecules e.g., based upon fit with 
another macromolecule - a receptor, enzyme, or in this case 
DNA, it is difficult and, in some cases, impossible to construct a 
reliable pharmacophore. This problem results in part from the 

25 fact that even closely related active molecules frequently fit into 

macromolecules in very different ways. 

A pharmacophore, as used herein, is defined as a 3- 
dimensional shape having a specific volume derived from the 
combined van der Waals surface of active molecules oriented by 
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fit into DNA, coupled with point charges located adjacent to the 
surface. A pharmacophore represents an aggregate array of 
positions in space of a series of molecules having die same or 
similar biological activity. The van der Waals surface can be 

5 represented in various ways including as a volume map, a dot 

surface, or a Connolly surface. The point charges are represented 
as dummy atoms whose positions are determined by die average 
positions of functional groups on active molecules which can form 
hydrogen bonds. Suitable charges are placed on die dunraiy atoms 

10 consistent with the capacity of the active molecules to form 

hydrogen bonds. The pharmacophores are specific for different 
compounds, their related molecules and a particular biological 
activity. According to the present invention, within the general 
class of molecules called hormones, an estrogen pharmacophore, 

15 an anti-estrogen pharmacophore, an androgen pharmacophore, a 

thyroid hormone pharmacophore, and a toxicity pharmacophore 
(shown in Figures 2-5 and Appendix 1) have been disclosed. It 
should be emphasized that these created pharmacophores do not 
exist as such in nature and are the product of aligning several 

20 related molecules to common binding sites in DNA using methods 

as described herein. 

Many other pharmacophores have been constructed 
using the method described in this application. A pharmacophore, 
once created, stands alone and is subsequently independent from 

25 the nucleic acid that was involved in its formation. Thus, after 

formation of the pharmacophore, one no longer needs to use the 
DNA as a template for the design of biologically active molecules. 
The pharmacophore itself can be used to generate new molecules 
that will possess the same or similar structural and charge features 
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that are represented by the pharmacophore. This is a completely 
different concept from the one of using the DNA as the model for 
the design of compounds. The pharmacophore can be used itself 
for any number of applications, including but not limited to the 
5 following: as a screening tool for drug development; to determine 

if a particular compound will possess bioactivity of a certain type, 
for instance estrogenic or androgenic activity; for toxicological 
evaluation; and to design compounds that possess increased or 
decreased binding affinity for DNA. 

10 Each pharmacophore has a characteristic shape, 

topology, volume^ and electrostatic profile. A pharmacophore is 
accurately described by its three dimensional shape which is 
represented by a coordinate system that is configured in computer 
memory (see Figures 2-5 for examples of pharaiacophores which 

15 are also contained in data files in Appendix 1 , the magnetic QIC 

tape submitted with the application). Each specific atom within a 
molecule that fits in a pharmacophore has a specific location 
relative to the docking heteroatoms. The individual atoms also 
have electrical charges assigned to them. These charges are 

20 represented numerically and through many other ways including 

the use of colors and shading to indicate field strength. As the 
degree of steric and electrostatic fit between the pharmacophore 
and the dummy atoms increases, resulting in a negative energy of 
interaction (-kcal), the efficacy of the pharmacophore increases 

25 which could manifest as increased bioactivity. The term "energy 

of interaction" as used herein is the total energy in -kcal of a 
molecule as it is being fitted into a pharaiacophore. This has been 
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observed in the case of molecules that fit within the estrogen 
pharmacophore and bioactivity in a uterotropic assay. The 
volume of a pharmacophore is described in cubic angstroms. The 
pharmacophore can be cross sectioned precisely in any plane and 
internal distances measured with an angstrom ruler. The 
circumference of any cross section is easily measured with 
morphometric analysis. Similarly, specific subregions of the 
pharmacophore, such as the site that binds to the DNA, can be 
subjected to the same methods of analysis. 



Construction and Utility of Pharmacophores 

An example of the construction and utilization of the 
estrogen pharmacophore is given below. Figure 2A is a computer 
generated space filling stereo view of the DNA cavity which fits 

15 estrogens. The fit of active estrogens oriented by energy 

calculations into the DNA cavity, using the methods described 
above, is presented in Figure 2B. Figiure 2C shows the combined 
active surface of estrogens removed from the cavity in DNA that 
is used to constmct the pharmacophore. The atoms are colored in 

20 the following manner: carbon/white; hydrogen/cyan; 

nitrogen/blue; oxygen/red; phosphorus/yellow. Figure 3 
demonstrates a volume contour map (yellow) in stereo with 
dunmiy atoms (magenta) surrounding the active molecules which 
were used in the construction of the pharmacophore (A); the 

25 empty pharmacophore (B); fit of the highly active estrogen 

3,1 lp,17p-Trihydroxy-7a-metiiylestra-l,3,5(10)-triene 1 1-nitrate 
ester (hereinafter 7a-methylestradiol-l Ip-nitrate ester reported 
in Peters et al., J. Med. Chem, 32:2306-2310 (1989)) which is 
acconunodated completely within the pharmacophore (C); poor fit 
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of the inactive estrogen 9|3-estradiol which extends appreciably 
beyond the surface of the pharmacophore (D). Figures 4 and 5 
present examples of the three dimensional appearance of estrogen, 
androgen, thyroid, anti-estrogen, and toxicity pharmacophores. 
5 The data files for the volume maps and dummy atoms for each of 

these pharmacophores and for Figures 2 - 5 are contained in 
Appendix 1, the magnetic QIC 120 tape submitted widi the 
application. This QIC computer tape has the following format and 
files generated from Sybyl 6.0 mn on a Silicon Graphics Indigo 
10 Extreme. 

1 ) /usr/people/guest/empty .dir/estrogempharmacophore2.dsp 
This is the volume map for the estrogen pharmacophore. 

2) /usr/people/guest/empty.dir/androgenpharmacophore.dsp 
This is the volume map for the androgen pharmacophore. 

15 3) /usr/people/guest/empty.dir/thyroid pharmacophore.dsp 

This is the volume map for the thyroid pharmacophore. 

4) /usr/people/guest/empty.dir/antiestrogenpharmacophore.dsp 
This is the volume map for the antiestrogen pharniacophore. 

5) /usr/people/guest/empty.dir/toxicitypharmacophore.dsp 
20 This is the volume map for the toxicity pharmacophore. 

6) /usr/people/guest/empty.dir/dummyatoms.mdb 

This is the molecular database file containing the appropriate 
dummy atoms for each pharmacophore designated by name. 

7) /usr/people/guest/empty.dir/patent94final 

25 This is a file containing the four color prints (Figures 2 through 

5) submitted with the patent application describing the 
pharmacophores constructed with the Silicon Graphics program 
Showcase. Copies of the Appendix 1 tape containing the data files 
are also retained by the inventor. These files contain the 
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information that represents the three dimensional configuration of 
each pharmacophore. These files may preferably be read with 
Sybyl software (version 6.03; Tripos Associates, St. Louis, MO) 
running on a computer system such as; CONVEX, ESV, IBM 
5 RISC/6000, IRIS 4D, SUN 4, or VAXA^MS, equipped with a QIC 

tape drive. In addition, other media such as 4 and 8 mm DAT 
tapes may be used to store and read these data files on the systems 
mentioned above. 

Quantitative measurements of the degree of fit of 

Fit is determined by measuring the amount of volume of each 
structure which could be placed within the pharmacophore 
volume map and normalizing the value to that of the natural 
hormone estradiol set at 50%. Electrostatic interactions with the 

15 dummy atoms are optimized for each compoimd and calculated 

using the Tripos force field. The electrostatic energy value is 
normalized to that of estradiol set at 50%. In this study, the total 
of volume fit and electrostatic fit are treated equally and totaled to 
reflect the overall fit in xhc pharmacophore. As shown in Figure 

20 6, degree of fit to the pharmacophore correlates highly to relative 

uterotropic (estrogenic) activity. In comparison to estradiol, 7a- 
methylestradiol-lip-nitrate ester (labeled 2), which is not part of 
the data set used to construct tiie pharmacophore, fits appreciably 
better than estradiol (labeled 1). In contrast, 7a-methylestradiol- 

25 11 p-nitrate ester binds poorly to the estrogen receptor (less than 

6% of the binding of estradiol). The uterotropic v£dues for Ta- 
me thy lestradiol- 11 P-nitrate ester (labeled 2) relative to estradiol 
set to a normalized value of 100% are considerably greater than 
that of estradiol (Peters et al., /. Med. Chem. 32:2306-2310 
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(1989)) as predicted by fit into the pharmacophore. In contrast, it 
is not possible to fit 9|3-estradiol (labeled 3) into the 
pharmacophore and, as predicted, this analog has little uterotropic 
activity. In summary, fit of compounds to the pharmacophore 
5 correlates with biological activity and can thus be used to design 

new compounds by virtue of their fit. It is noteworthy that the 
highly potent estrogenic nitrate ester binds very poorly to the 
estrogen receptor. Thus, it would not be possible to predict the 
estrogenic activity of this analog on the basis of receptor binding 

a protein receptor. In many cases, compounds with greater 
estrogenic activity than the natural hormone estradiol bind 
relatively poorly to the estrogen receptor. 

The fit of compounds into DNA is consistent with, 

15 but not the same as, fit into the pharmacophores. The distinction 

between the DNA cavities and the pharmacophores is that the 
surfaces were derived from different structures, i.e. the DNA 
cavities from the DNA structure and the phamaacophores from 
the combined surfaces of active compounds. As shown in the 

20 examples and in Figures 2-5, and in datafiles contained on the 

magnetic tape as Appendix 1, the degree of fit to the 
pharmacophores is obtained by fit to the three dimensional map 
which represents the active compounds. The pharmacophores 
enable quantitative determination of the degree of fit to the 

25 combined surfaces of the active compounds and this information 

cannot be obtained from fit into DNA. Moreover, the fit of 
active compounds to the pharmacophore can be quantitated based 
upon the portion of the molecule which does not fit within the 
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pharmacophore volume. This enables automatic assessment of 
predicted inactive structures. 

Toxicity Pharmacophores (Toxicophores) 

Pharmacophores can be constructed to represent a 
three dimensional shape that is predictive of toxic biological 
activity. Such pharmacophores, called toxicophores, have regions 
that would potentially damage DNA. A typical toxicophore has 
been constructed using, tetrodotoxin, dioxin, RU 486, dilantin, 
thalidoniide and oroflex, among other compounds. An exam.ple of 
this toxicophore is provided in Figures 4E and 5E and on the 
submitted magnetic data tape as Appendix 1. For example, by 
overlaying this toxicophore on another pharmacophore such as 
the estrogen pharmacophore (Figure 4A and 5A), a drug designer 
would know to avoid designing an estrogenic compound with 
certain molecular groups that might impart toxic activity if these 
groups extended into the three dimensional space occupied by the 
toxicophore. This approach would greatly facilitate and 
economize drug design by guiding the designer to avoid 
synthesizing estrogenic compounds that might have damaging 
effects on DNA as opposed to proceeding with synthesis and 
purification and subsequently discovering that the compound 
possesses dangerous toxicity. 

Solvent Pharmacophores (Aquaphores) 

Pharmacophores in their relationship to nucleic acids 
are usually surrounded by a solvent. The predominant solvent in 
living organisms is water and accordingly, most pharmacophores 
exist in an aqueous environment. Water is the preferred 
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embodiment of the solvent pharmacophore and is termed an 
aquaphore. Pharmacophores, and their molecules may also be 
placed in non-aqueous environments for various purposes such as 
crystallographic studies or other analytical procedures. 
5 In living organisms, the aqueous environment 

surrounding the pharmacophore also has an intimate association 
with the adjacent nucleic acid. This aqueous shell assists in the 
optimal fit of the pharmacophore into the cavity of the double 
stranded DNA, and has its own three dimensional shape. The 

10 optimal steric and electrostatic placement of water molecules in 

tiie space between the pharmacophore and the DNA is achieved in 
the present invention. This three dimensional shape is called a 
solvent pharmacophore, and can be described in all the ways Usted 
above for the pharmacophores based on oUier molecules such as 

15 estrogen. Solvent pharmacophores assist the designer of 

compounds by placing limits on the dimensions of a compound 
designed using a particular pharmacophore as a template. In 
addition, solvent pharmacophores assist the creator of 
pharmacophores because the solvent shell or cage represented by 

20 the pharmacophore provides enhanced ability to properly align 

molecules relative to DNA during the creation of the 
pharmacophore. 

Receptophore 

25 Many molecules, such as steroid hormones, are 

shuttled to the nucleus by other molecules known as receptors 
(Tsai and O'Malley, Ann, Rev, Biochem. 63:451-486 (1994)). 
These receptors bind tiie hormones (called Ugands), bind to the 
nucleic acids, for example in their DNA binding domain, and 
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present ligands to nucleic acids such as DNA. Evidence suggests 
that the binding of the receptor to the DNA causes a 
conformational change in the DNA to facilitate insertion of the 
ligand (NarduUi et al., Molec. Endocr, 7:331-340 (1993)). The 
5 pharmacophore concept is based on the three dimensional shape of 

the optimal fit of related molecules into nucleic acids such as 
partially unwound, double-stranded DNA. The DNA binding 
domain of the receptor can be modeled into a three dimensional 
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Example 1 

Fit into partially unwound double stranded DNA 
using the mammalian female hormone estradiol. 
5 Computer modeling has demonstrated that the 

mammalian steroid progesterone is a remarkable "lock and key" 
fit into DNA at 5'-dTdG-3', 5'-dCdA-3' (Figure 1). Each of the 
known x-ray crystal structures of progesterone is capable of 
forming two stereospecific hydrogen bonds and a stable complex 

10 measured by force-field calculations. Remarkable 

complementarity is evident in the complex by the overlap of 
hydrophilic and hydrophobic regions of the steroid and DNA. 
The enantiomer of progesterone which does not occur in nature 
does not fit. The plant hormone gibberellic acid has also been 

15 shown to fit between base pairs but in a different sequence, i.e., 

5'-TdA-3', 5'-dTdA-3'. Four stereospecific hydrogen bonds are 
formed within the couples: ent- gibberellic acid does not fit. 

The mammalian female hormone estradiol also fits in 
DNA (Figure 1). Two stereospecific hydrogen bonds of 

20 approximately 2.65 A are formed between each hydroxy 1 group 

of the steroid and phosphate oxygens on adjacent strands. The 
overall fit within the complex is about -59 kcal. Most alterations 
of the positions of the hydroxyl groups on the estratriene nucleus 
result in a substantial loss of potential electrostatic interactions 

25 with the DNA. Moreover, most alterations of the absolute 

stereochemistry of the cyclopentanophenanthrene ring pattern also 
result in a substantial loss of potential, electrostatic interactions 
with the DNA. Moreover, most alterations of the absolute 
stereochemistry of the cyclopentanophenanthrene ring pattern also 



SUBSTITUTE SHEET (RULE 26) 



wo 95/14791 



PCTAJS94/13765 



38 



result in a poor fitting molecule. This is evident in attempting to 
fit 9P-estradiol into DNA. The puckering in the steroid caused by 
inversion of the estradiol stereochemistry at C-9 from a to (3 
prevents complete insertion between the base pairs. Even if 
5 forced into DNA without regard to strain caused by the overlap 

of van der Waals surfaces, 9(i-estradiol can form only a single 
hydrogen bond. The relative fit into DNA resulting from partial 
insertion of 9p-estradiol is about -17 kcaL 9p-estradiol is inactive 
when tested in vivo for estrogenic (uterotropic) activity. 
10 The finding that estradiol is a "lock and key" fit into 

DNA, although most structural alternatives to estradiol fit poorly 
provides further support for the premise that DNA 
stereochemistry contains the master blueprint for natural product 
structures. 

Example 2 

Correlation of Estrogenic Activity with the Fit of 
Estrogens and Related Analogs into DNA. 

That fit into DNA measured by energy calculations 

20 can be correlated with biologic activity was demonstrated using a 

series of estrogens and related synthetic analogs. The molecules 
which are inactive in uterotropic assays fit poorly into DNA. 
Molecules that fit into DNA better than estradiol are all more 
active than estradiol in the uterotropic assays, for example, lip- 

25 acetoxyestradiol (approximately - 68 kcal). This correlation is 

also observed with nonsteroidal, synthetic estrogens, such as the 
potent synthetic estrogen, transdiethylstilbestrol (approximately - 
62 kcal), which fits well, whereas the poorly active geometrical 
isomer cis-diethyistilbestrol (approximately -20 kcal) is a poor fit. 
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Example 3 

Biosynthetic Pathways Reflect Increasing Fit into 
DNA Whereas Inactivation Pathways Lead to Decreasing Fit into 
5 DNA. 

Molecular modeling studies conducted with the 
mammalian hormone progesterone and the plant hormone 
gibberellic acid have shown that each step in the respective 
biosynthetic pathway reflects a structural change that results in 

10 increased fit in DNA. For example, when considering possible 

stereoisomers which could result in any given step in 
progesterone biosynthesis, the best fitting structure is one which 
was produced in nature. In sharp contrast, each step in the 
inactivation of progesterone eventually leading to the excreted 

15 gluciuronides and sulfates resulted in the worst possible fitting 

stereoisomers. 

The two possible dihdro reduced metaboUtes of the 
male hormone testosterone were examined for fit and correlation 
with biological activity. Relative to testosterone (100%), 5p- 

20 dihydrotestosterone is a poor fit (84%), whereas its epimer 5a- 

dihydrotestosterone fits even better than testosterone (102%). 
These data are consistent with published findings by Hilgar and 
Hummel, "The androgenic and myogenic evaluation of steroids 
and other compounds -assay 1", in AG Hilgard. D.J. Hiuiunel 

25 (ed.) Endocrine Bioassav Data. Part III )U.S. Dept. HEW NIH 

1964), that 5a-dihydrotestosterone is a highly active androgen, 
whereas 53-dihydrotestosterone is essentially inactive. 

In the case of the tetrahydro reduced stereoisomers 
of progesterone, 3a-hydroxy-5P-pregan-20-one was a poor fit 
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whereas its stereoisomer 3a-hyciroxy-5a-pregnan-20-one was an 
excellent fit. The former molecule is a highly active neurosteroid, 
as reported by Purdy, et al., /. Med, Chem. 33:1572-1581 (1990). 
Comparison of the specific pattern of donor/acceptor linkages of 
3a-hydroxy-5a-pregan-20-one with those in the steroid/thyroid 
hormone superfamily demonstrated that the linkage pattern of 3 a- 
hydroxy-5a-pregan-20-one is unique. Compounds having such 
unique linkages are predicted to have unique biologic function and 
may be ligands for newly discovered "orphan receptors/' 



Example 4 

Correlation of Toxicity and "Side Effects" with Fit 

into DNA, 

Ligands that fit into more than one site in DNA have 
15 been observed to have multiple biologic actions. Both desirable 

and undesirable "side effects" should thus be predictable from the 
specific DNA sequence which a given compound fits into as well 
as the manner and relative degree of fit. Examples of molecules 
that have been observed to fit into more than one site in DNA 
20 include the psychotropics cocaine, morphine, LSD and 

tetrahydrocannabinoids, and certain intercalating antibiotics. For 
example, the monoamine oxidase inhibitor selegiline fits into the 
site in DNA which accommodates glucose and various oral 
antidiabetic dmgs as reported by Rowland et al., 7. Clin, 
25 Pharmacol. 34:80-85(1994). This observation is consistent with 

the finding that selegiline causes hypoglycemia in some patients. 
Another case is the anti-androgen anandron which fits into the site 
in DNA which accommodates testosterone. Because anandron fits 
into DNA in two orientations, i.e., in a manner similar to both 
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androgens and anti-androgens, mixed activity is predicted for this 
compound. Experimental results indicate that anandron has both 
agonist and antagonist activities, as reported by Steinsapir, et aL, 
The Endocrine Society (74th Annual Meeting) 1992:109 (abs. 
5 228). 

It has also been noted that Ugands which cause stress, 
chemical modifications and/or covalent linkages to the DNA when 
fit into a given site frequently possess toxicity. Examples include 
certain carcinogens and teratogens, e.g., thalidomide, dioxin, 

10 arene oxides, aflatoxins and some diethylstilbestrol metabolites. 

Anotiier example is the anti-progestin RU486 which stresses base 
pair hydrogen bonds when inserted fully into the progesterone 
site in DNA. Similar strain is produced by other anti-progestins 
having the same side chains (e.g., 1 Ip-phenylamines) raising the 

15 possibiUty that such features may correlate with abortifacient 

activity attributed to RU486 and related analogs. Observations 
with thalidomide enantiomers indicate that teratogenicity 
associated with this compound may correlate with a stereospecific 
effect on base pairing. An example of such a toxicity 

20 pharmacophore, called a toxicophore, is presented in Figures 4E 

and 5E and the corresponding data file is submitted on magnetic 
tape. 



Example 5 

25 Design and Development of New Drugs. 

Using the principles described above and in the 
examples, new drugs can be designed or existing drugs can be 
redesigned while at the same time limiting potential undesirable 
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side effects. One example of how an active drug can be designed 
using the technology follows; 

3-PhenylacetyIamino-2,6-piperidinedione (AlO) is a 
modified amino-acid derivative, which was originally isolated 
from freeze-dried human urine. Despite having low toxicity, 
high concentrations of AlO were required to demonstrate 
significant growth inhibitory activity on tumor cells. The focus 
of the following study was to develop more potent analogs. 
Modeling studies demonstrated that AlO was capable of inserting 
into partially unwound double stranded DNA and forming a 
single hydrogen bond between the imino proton of the 
piperidinedione ring and a phosphate oxygen on a single strand. 
It was observed that placing a hydroxyl group at the para position 
of the phenyl group of AlO would enable a second hydrogen bond 
to form thereby substantially enhancing fit as reported by Hendry 
et al., y. Steroid Biochem. Molec. Biol 48:495-505 (1994), die 
teachings of which are hereby incorporated by reference in their 
entirety. The relative fit of AlO (normalized to 100%) and 
various related analogs measured by energy calculations 
demonstrate that the best fitting compound is p-OH-AlO (164%). 
Subsequent synthesis of these compounds followed by testing in 
various animal and human tumor cells demonstrated that p-OH- 
AlO was the most active compound and was as much as an order 
of magnitude more active than AlO, as reported by Hendry, et al.. 
Recent Advances in Chemotherapy, Buchner and Rubinstein (eds) 
2498-2499 (1991), Hendry et al., U.S. Patent 5,238,947 which is 
incorporated herein by reference. 
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Example 6 

Comparison of the Drug Design Technology to 
Classical Structure-Activity Methods, 

The drug design technology described here can be 
5 used in conjunction with quantitative-structure-activity- 

relationship methods (QSAR), e.g., comparative field molecular 
analysis (CoMFA). One value of the approach is that it facilitates 
the orientadon of various ligands relative to one another in three 
dimensions. The successful structure-activity relationship found 
10 for estrogens derived from fit into DNA is described here. If one 

were to attempt to derive such a relationship a priori without first 
knowing the detailed three-dimensional structure of an 
appropriate macromolecule (e.g., the ligand binding site of a 
receptor or an enzymatic site), chemical intuition would 
15 necessitate searching for conmion features that exist in known 

active stmctures. In the case of the natural hormone estradiol and 
the potent synthetic estrogen trans-diethylstilbestrol, such a 
common feature is the phenoxy group. Alterations of the 
phenoxy group give rise to inactive structures. If one overlaps 
20 the resulting three dimensional orientation with that which is 

obtained by optimal docking of these molecules into DNA, a 
different pattern emerges. Thus, one would expect that such 
different orientations when used subsequently to correlate activity 
of other molecules would give very different results. In fact, 
25 analysis of relative fit into DNA in kcal shows that if the 

orientation of trans-diethylstilbestrol based upon overlapping the 
A ring of estradiol is used to dock die ligand into DNA, a 
relatively poor fit results. In this case, poor activity for trans- 
diethylstilbestrol would be incorrectly predicted. In contrast. 
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using the orientation of DES derived only from the 
stereochemistry of DNA as taught by the present invention, 
increased activity would be correctly predicted. 

Example 7 

Design of an anti-estrogen, para- 
hydroxyphenylacetylamino-2, 6'piperidinedione, a re gulator of 
tumor cell growth. 

Within the past few years there has been a growing 
interest in nontoxic, naturally occurring small molecules as 
regulators of tumor cell growth. Examples of recently published 
findings include: regression of mammary carcinomas by a dietary 
monocyclic nonoterpens, limonene; modulation of oncogene 
expression in erythroleukemic cells growth by an endogenous 
product of lipid peroxidation, 4-hydroxynonenal; inhibition of 
malignant cell growth by the endogenous ligand p- 
hydroxyphenylacetate; and induction of tumor cell differentiation 
in premalignant and malignant cells by a circulating component of 
human plasma, phenylacetate. 

Phenylacetate has been shown to reduce levels of the 
myc oncogene which is involved in the development of several 
cancers including breast, brain, prostate, blood, lung and colon. 
Another mechanism by which phenylacetate is thought to be 
effective is by reducing levels of the amino acid glutamine. 
Phenylacetate conjugates with circulating glutamine to produce 
the excreted urinary metabolite phenylacetylglutamine (PAG). 
Cancer cells require glutamine for growth and are known to be 
more sensitive to glutamine depletion that normal cells. These 
findings have led to the initiation of Phase I clinical trials with 
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phenylacetate in brain and prostate cancer at the National Cancer 
Institute. 

In the process of screening fractions of freeze dried 
human urine for growdi inhibitory activity in human breast 
5 cancer cells, a dehydration product of PAG was isolated. The 

compound was characterized as 3-phenylacetylamino-2,6- 
piperidinedione by spectroscopic methods and independent 
synthesis and was termed A 10 based upon the chromatographic 
fraction from which it was isolated. It has not been conclusively 

however, it is similar to phenylacetate in that it lacks toxicity in 
both laboratory animals and humans. The compound has also 
been found to inhibit growth in a variety of cancer cells in vitro, 
as well as human breast cancer transplanted into athymic nude 

15 mice. Chemoprevention effects have also been reported. In 

addition to phenylacetate and methyl p-hydroxyphenylacetate, 
there is a number of synthetic compounds with antitumor 
activities that have structural features in common with 3- 
phenylacetylamino-2,6-piperidinedione. Examples of such 

20 compounds in which analogies to the piperidinedione ring are 

prominent include the aromatase inhibitors aminoglutethimide, 
rogletimide and related analogs which are used in the treatment of 
breast cancer; the alkylating drug PCNU which inhibits tumor 
growth by proposed interaction with DNA; 5-cinnamoyl-6- 

25 aminouracil derivatives which inhibit tumor growth via putative 

DNA intercalation; amonafide and its congeners which mediate 
topoisomerase II DNA cleavage by intercalation, bis (2,6- 
dioxopiperazine) derivatives, e.g., ICRF-193, which are potent, 
direct inhibitors of mammalian DNA topoisomerase II. 
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Relative high doses of 3-phenylacetylainino-2,6- 
piperidinedione have been generally required to inhibit tumor 
growth both in vitro and in vivo. The goal of this study was to 
identify more active phenylacetylamino-2, 6-piperidinediones 
5 using the technique described above. The technology is based 

upon modeling of the stereospecific fit of molecules into DNA 
and has been recently modified to take advantage of computer 
graphics and energy calculations. Computer modeling was 
followed by the design, synthesis, and in vitro biological testing of 

10 3-phenylacetylamino-2,6-piperidinedione derivatives. The 

molecule predicated to be the most active based upon degree of fit 
in DNA, i.e., die p-hydroxy derivative, was found to be the most 
active antitumor agent in all of the biological assays investigated. 
When tested in MCF-7 (E3) human breast cancer cells, the p- 

15 hydroxy derivative possessed antiestrogenic activity in the range 

of the drug tamoxifen which is currently in clinical use for the 
treatment of breast cancer. 

Example 8 

Three separate molecules designed with modeling 
to technology are shown in Figure 7. These molecules, termed SGI 

100, SGI 101 and SGI 102 bear similarities and differences to 
each other. They all show structural similarities to components of 
both estrogen and progesterone. SGI 100 was designed on the 
basis of its ability to fit into the site in DNA which accommodates 
5 both estradiol and progesterone. The manner in which it fits 

predicted antagonist activity. When fit to the estrogen 
pharmacophore (Figure 2 and Appendix 1), die acetyl group at 
the np position extended out of the pharmacophore and had 
electrostatic repulsion with dummy atoms (-131 kcal) compared 
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to estrogen (-51 kcal). Accordingly, this high positive energy of 
interaction indicates that SGI 100 acts as an antagonist. The 
binding of SGI 100 to the estrogen receptor is dose dependent and 
approximately 144 times less than estradiol. In bioactivity 
5 experiments, SGI 100 significantly decreased cell growtii in MCF- 

7 human breast cancer cells (134,431 cells) at a dose of 10*^ M 
when compared to control cells (252,197 cells). The same 
concentration of tamoxifen citrate decreased the number of MCF- 
7 cells to 187,759. Thus, the design of this compound based on the 

10 pharmacophore approach of this invention predicted a 

demonstrable anti-estrogen bioactivity that was greater than 
tamoxifen. SGI 100 binds in a dose dependent manner to the 
progesterone receptor but with 133 to 200 times less affinity. 

SGI 101 was designed on the basis of fitting into 

15 DNA at the site which accommodates estradiol but with opposite 

hydrogen bonding properties which predict estrogen antagonist 
activity. SGI 101 extends beyond the estrogen pharmacophore and 
has electrostatic repulsion between the para-nitro group and 
dummy atoms. SGI 101 is the most potent analog designed by the 

20 technology as measured by growth inhibition of MCF-7 cells. At a 

dose of 10"** M, SGI 101 inhibited cell growth (81,103 cells) 
relative to control (252,197 cells) and was substantially more 
active than the same concentration of tamoxifen (187,759 cells). 

SGI 102 was designed on the basis of its fit into DNA 

25 at the site which acconimodates progesterone. SGI 102 possesses 

an alkyl amino side chain at the 1 ip position which extends out of 
the site between base pairs into the major groove. It has different 
hydrogen bonding properties than progesterone and would extend 
beyond the volume map of the progesterone . pharmacophore. As 
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such, it is predicted to be an antagonist. SGI 102 was designed 
prospectively, synthesized, and tested in various biological assays. 
SGI 102 binds in a dose dependent manner to the progesterone 
receptor but not as strongly as progesterone or the abortifacient 
5 antiprogestin RU486. In animal experiments, SGI 102 showed no 

abortifacient activity. However, in experiments using MCF-7 
human breast cancer cells, SGI 102 had equivalent activity to 
RU486 in inhibiting growth. These findings are consistent with 
the predictions made by the modeling technology. 
10 It should be understood that the foregoing relates 

only to a preferred embodiment of the present invention and that 
numerous modifications or alterations may be made therein 
without departing from the spirit and the scope of the invention as 
set forth in the appended claims. 
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CLAIMS 

I claim: 

1. A method for designing molecules having 
5 biological activity comprising 

designing compounds which are complementary to a 

double stranded DNA helical structure; 
determining the energy of interaction between the 
compound and the DNA in the presence of specific 

selecting compounds having an equal or lower energy 
of interaction than the natural ligands for the 
specific biological receptors, 

^5 2. The method of Claim 1 wherein the receptor 

binds to a compound selected from the group consisting of 
steroids, thyroid hormones, vitamins, phytohormones, peptide 
hormones, neurotransmitters, prostaglandins, and sugars. 

20 3. The mediod of Claim 1 wherein the compound 

has a cyclopentanophenanthrene motif. 

4. The method of Claim 1 wherein the compound 
is complementary to an unwound cavity within a DNA helix. 

25 

5. The method of Claim 1 wherein the compound 
is complementary to a cavity in DNA formed by shifting a base in - 
an unwound site. 
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6. The method of Claim 1 wherein the compound 
is complementary to a cavity formed in DNA by removing a 
nucleotide base. 

5 

7. The method of Claim 1 wherein the compound 
is a receptor antagonist. 

8. The method of Claim 7 wherein the compound 
10 is reproductive hormone antagonist. 

9. A compound having biological activity isolated 
by the process comprising, 

designing compounds which are complementary to a 
15 double stranded DNA helical structure; 

determining the energy of interaction between the 
compound and the DNA in the presence of specific 
biological receptors; and 
selecting those compounds having an equal or lower 
20 energy of interaction than the natural ligands for the 

specific biological receptors. 

1 0. The compound of Claim 9 wherein the 
receptor binds to a compound selected from the group consisting 

25 of steroids, thyroid hormones, vitamins, phytohormones, peptide 

hormones, neurotransmitters, prostaglandins, and sugars. 

11. The compound of Claim 9 wherein the 
compound has a cyclopentanophenanthrene motif. 
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12, The compound of Claim 9 wherein the 
compound is complementary to an unwound cavity within a DNA 
helix. 

13. The compound of Claim 9 wherein the 
compound is complementary to a cavity in DNA formed by 
shifting a base in an unwound site. 



compound is complementary to a cavity formed in DNA by 
removing a base. 



15. The compound of Claim 9 having estrogenic 



15 activity. 



16. The compound of Claim 9 wherein the 
compound is a receptor antagonist. 

20 17. The compound of Claim 16 wherein the 

compound is a reproductive hormone antagonist. 

18. A pharmacophore comprising, 
a three-dimensional array of points defining a specific 
25 shape £uid volume, wherein said three-dimensional 

array of points is the aggregate average shape of a 
plurality of molecules when optimally fit into an 
unwound DNA site, said molecules possessing the 
same or similar biological activity. 
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19. A biologically active compound designed using 
a pharmacophore of Claim 18, wherein the template identifies 
areas of structure and charge responsible for fit into DNA. 

20. A method for creating a pharmacophore 
comprising the steps of: 

determining the optimal fit of a plurality of 
compounds having the same or similar biological 
activity into nucleic acid sequences such that the 
lowest energy of interaction and best steric fit are 
obtained; 

aligning the compounds relative to the heteroatoms on 

the nucleic acids; and 
defining a three dimensional shape representing the 

aggregate average shape resulting from said 

alignment of the compounds, 

21. The method of claim 20, wherein the nucleic 
acid sequence is selected from the group consisting of 
deoxyribonucleic acid, double stranded deoxyribonucleic acid, 
ribonucleic acid, complexes of deoxyribonucleic acid and 
ribonucleic acid, and apurinic and apyrimidinic sites. 

22. The method of claim 20, wherein the nucleic 
acid is deoxyribonucleic acid. 
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23. A method of screening a molecule for 
bioactivity or toxicity wherein the molecule is examined for the 
ability to fit into a pharmacophore comprising the steps of: 

measuring the energy of interaction when the molecule 
5 is fitted into the pharmacophore; and 

comparing the energy of interaction of the molecule to 
a predetermined energy of activation that correlates 
to a biological activity. 

pharmacophore is selected from the group consisting of an 
estrogen pharmacophore, an antiestrogen pharmacophore, an 
androgen pharmacophore, a thyroid hormone pharmacophore, 
and toxicophore. 

15 

25. A method for screening a molecule for toxicity 
wherein the molecule is examined for the ability to fit into a 
toxicity pharmacophore comprising the steps of 

measuring the energy of interation when the molecule 
20 is fitted into the pharmacophore; and 

comparing the energy of interaction of the molecule to 
a predetermined energy of activation that correlates 
to a toxic activity. 
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26. A method of designing a molecule with a 
desired biological activity comprising the step of determining the 
stmcture of a molecule which will have a minimal predetermined 
energy of activation when fitted into a pharmacophore with a 
5 desired activity. 
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